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SECTION I 
INTRODUCTION 


During a Shuttle launch, structures in the proximity of the launch pad are 
subjected to acoustic pressures generated by rocket exhausts. The design of 
some structures, particularly those having a large area-to-mass ratio, is 
governed by the launch-generated acoustic environment, a relatively short but 
very intense pressure transient. 

A procedure documented in KSC-DM-3147, Procedure and Criteria for Conducting a 
Dynamic Response Analysis of Orbiter Weather Protection System on LC-39B Fixed 
Service Structure, was developed to calculate the peak dynamic response of a 
normal vibration mode of a structure to an input transient. The procedure is 
based on modal parameters of the structure, response spectra to acoustic pres- 
sures, and a definition of acoustic field by means of pressure correlation 
lengths. Modal parameters and the area of correlated pressures bounded by 
pressure correlation lengths are a prerequisite to calculations of generalized 
modal loads required as input to response spectra for a consequent estimate of 
peak responses. 

Generalized modal loads are defined by means of AJ-coefficients, which are a 
measure of a vibroacoustic coupling between the structure and the acoustic 
field. The stronger the vibroacoustic coupling is, the higher the vibration 
of the structure in the coupled mode, and vice versa. Thus, an accurate esti- 
mate of AJ for each vibration mode is essential to calculations of dynamic 
response. Even for a structure of an average complexity, calculations of AJ- 
coefficients may become a cumbersome task because existing dynamic analysis 
codes do not have the capability to automatically account for all parameters 
leading to the computation of the required extreme values of AJ's. These 
coefficients depend on a surface integral of the product between modal dis- 
placements and the correlated pressure distribution, a function of pressure 
correlation lengths that depend upon a particular acoustic field and the 
resonance frequency of the mode. The relative position between the mode and 
the center of distribution affects variable integral limits and the values of 
AJ-coefficients. 

A position resulting in the extreme value of AJ and in the highest dynamic 
response of a mode mey not be obvious from an examination of integrand con- 
stituents and variable integral limits. Particularly in vibration modes where 
modal displacements of the surface under an acoustic load change signs in re- 
lation to the outward surface normal, a proper optimization of the load posi- 
tion becomes critical and, if not performed, may lead to significant dif- 
ferences (by factors exceeding 100 percent) in the estimate of a peak dynamic 
response. Generally, a different position of correlated pressure distribution 
corresponds to each vibration mode, resulting in the strongest vibroacoustic 
coupling and the highest response. 
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An example of the required op fi rrn zati on ^g^and/o^rcas^of the third 
fundamental mode of a three-span conti n b / when the pressure 

mode of a simply supported beam illustrated in rigu ^ ^ ^ f . ye 

correlation length in the ° h the P position of the distribution 

times the distance between note} the estimate of AJ-coeffi- 

center relative to the mode becomes c ' t 1 t center of correlated pressure 

cients. The extreme value of ^ near the end nodal point (the end sup- 

distribution (load center) is loc certainly not apparent, 

port of a three-span beam), a load posit . , center j s located anywhere 

If the load position is not optimized and and wit hin the 

near and between the second and t va l ue 0 f AJ becomes substantially lower 

than* the^extreme ia^'lenerally’belcjj^O Pf^^^^^Sts may'rfs^t 

|«^ u ^Sse p desi 9 n safety 
factors are usual ly Saar or below a factor of two. 

This example of h W sensitive compotatiohs^f 
solution and the output of a dynamic analysis code. 

This report presents diagrams for sl^le^d jltispan conti nuous^eams^from 

v^lid application -e stated i n t e presen ed f de^i vat, on. ,^1 the^exp^ 
by means of cited assumptions requ v frpnupncv band acoustic 

¥?;’, e d "stress-strain ■ -"I* occ2r mainly in a single 

fundamental mode of each structural component. 

Applications may be extended to a grid of beam-type structures and t< . plates. 

-ductUVf 

“ ff a e r * ?? .YS'X't 2* or 1 n t fr^oriion ti t^ number^f .continuous 
“al^Thui? a mSl^fJfd reduction of response can be achieved by imple- 
menting a structural continuity. 
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The beam Is assumed to have a uniformly distributed mass (m) and constant 
stiffness (El). 



SYMMETRY (Load positions are 
from t symmetry only.) 


shown on the left 



First mode of a three-span beam 
(1 * span), or a 

third mode of a simply supported 
beam with a span L = 31 . 



Pressure 


31 < X<51 

Correlation Length 


Position of correlated pressure dlst 
ribution (CPD) yielding an extreme 
AJ -coefficient 


CPD CENTERS 



Positions of correlated pressure 
distributions (0.881<e <1.51). 
Yield AJ-coefficients which are 
substantially lower than the 
extreme. Differences may exceed 





31 < X< 51 


J 


Horizontal scale represents a case X s 41 


Figure 1-1. 


Effect of Load Position (CPD Center) on AJ-Coefficients, 
a Measure of Vibroacoustic Coupling 


1-3/1— 4 



KSC-DM-3265 


SECTION II 
PROBLEM DEFINITION 


Utilization of response spectra, Y=q/(AJ/M/u; n ), to acoustic pressures, p(t), 
in an application to the analysis of peak structural responses, requires a 
computation of a generalized modal load, G(t)=AJp(t), for each normal mode of 
a vibrating structure. The time history of acoustic pressures, p(t), is 
assumed to be known from measurements taken in the acoustic field where the 
structure is located. Generally, a multitude of measurements is required for 
a proper definition of all basic parameters of an acoustic field. In the 
following equations, it is assumed that measurements are available and 
required parameters are defined, so that the mathematical formulation can 
reference them as if they were obtained from a single set of measurements. 

Response spectra and pressure correlation lengths (PCL's) are assumed to be 
available for the frequency range containing resonance frequencies of all 
normal modes considered in the analysis. Computations of generalized modal 
loads are then reduced to the problem of estimating AJ-coefficients for each 
normal mode, and peak response modal coordinates, q's, are calculated from 
response spectra: 

2 

q = Y(AJ/M/cj n ) 

The equation defining AJ for a mode of a planar structure in x,y plane is: 

AJ = f 4»(x,y)P(x,y)dA (1) 

Ai 

where: 

4»(x,y) is normal mode displacements perpendicular to the x,y plane. 

P(x,y) is the general form of the distribution of correlated pressures 

in the x,y plane and within the correlation area, Ap. Implic- 
itly, distribution is also a parametric function of the reso- 
nance frequency, fn, of the mode <t>(x,y). Condition 0<P(x,y)<l 
limits the values of the distribution. 

Ai is an area defining the extent of the integral. It is equal to 

the overlapping area between the pressure correlation area, Ap, 
and the structural area. As, exposed to acoustic pressures. By 
defining an overlap, Ai strongly depends on the position of the 
center of correlated pressure distribution relative to the 
structure and, thus, relative to modal displacements. 

The coordinates of the distribution center, a and b, satisfy condition 
P(a,b)=l. However, a and b are also independent parameters that must be 
varied to yield an extreme value of AJ. Therefore, P(x,y) must be defined 

2-1 
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relative to its center while, in a homogeneous acoustic field, the center may 
£.™luio£ed anywhere on the surface of the structure. In the local coordi- 
nates u=x-aandv“y-E! the definition of distribution relative to Us center 
becomes: 


P(x,y) = P(a+u,b+v) = P uv (u,v) 
with limits on variation of u and v: 

-A x /2 < u < A x /2 and -X y /2 < v < X y /2 


( 2 ) 


(3) 


where X x and X y a™ pressure correlation lengths in the direction of the x and 


(4) 


y axes. 

The equation: 

P uv (u,v) = 0 

defines a contour of the pressure correlation area, Ap, in tjf x ,y plane. 
Distribution P, JV (u,v) is zero everywhere on and outside of this contour. 
Because pressure correlation lengths, A x and X y , 

from measurements of coherence and phase between available se "^ p * ;* ee 
KSC-DM-3147 Appendix B) and because the number of measurements is limited, 
«HflcUio; oflhe coniour shape where P uv (u,v)-0 was not poss,b •• 
tnur is assumed to be an ellipse with axes A x and X y . The elliptic contou 
introduces a constraint on the allowable variation of local coordinates u an 

v that is more stringent than (3): 

2 2 2 2 2 2 


4XyU + 4A x v < X x X y 


(5) 


Consequently for any point with coordinates u,v, new variables (q, X, r) are 
defi ned: 


q = ^ = X x X y 

■_'y ' 

Note that: for v = 0, X = A x 


1 + q‘ 


2 - q 2 X? 


1/2 : r = ±(u 2 + v 2 ) 1 / 2 


( 6 ) 


for u = 0, lim X-X y 


u -* 0 


Then, correlated pressure distribution (CPD) is explicitly defined. 
P uv (u,v) = \ + cos ^T L ) > 0 


(7) 


with a constraint 

- - < r < A 
2 2 


( 8 ) 
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Which ensures 
pressure area 


that point u 
Constraint 


,v remains inside the contour of the correlated 
(8) is equivalent to and replaces (5). The sign 


of 


r in 


(6) is irrelevant to the definition of P U v(u,v) 


The contour of the correlation area was assumed to be an ellipse resulting in 
the smaHestvalue of the correlation area (min Ap), but not necessarily , a 
lowering of AJ in (1). An assumption of a rectangular contour with sides x 
andX y results in the largest value of the correlation area (max Ap). The 
corresponding distribution becomes: 


Puv(^»^) 



+ cos 



+ cos 


2irv \ 

*y/ 


(9) 


with the constraint on the variation of u and v provided by (3). The form o 
(9) is equivalent to a definition of the CPD by means of two separate distri- 
butions (separation of variables) along each coordinate axis: 


P uv (u,v) - P u (u)P v (v) 


( 10 ) 


Bracketing of the extents of CPD's by (7) and (9) is not expected to have a 
significant effect on AJ. 

At this point, all constituents of AJ in (1) are formally defined, ^(x^y) Is 
known from either the finite element model (FEM) or from experimental modal 
analysis?" VulTuW) is defined by either (7) or (9) and by corresponding con- 
straints in (8) or (3). By locating the center of CPD at a point S(a,b) and 
by applying a coordinate transformation 

K (ID 

u = x-a , v = y-b 

to either 4>(x,y) ^ uv (u,v) or to P uv (u,v) -P(x,y ) so that both functions and 
their constraints are expressed in the same coordinates, the extent of the 
integral, the area Ai , can be defined. An additional constraint is imposed on 
the allowable variation of integrand coordinates by the choice of the point 
S(a b). This constraint contains previously defined (8) or (3) by requiring 
that integrand coordinates be contained within the contour of the overlap area 
Ai. An illustration of this problem is shown in figure 2-1. 

The integration of (1) can be performed in either the x,y or u,v coordinates 
and a value of AJ can be calculated. This value becomes a function of coor- 
dinates of the point S(a,b) that was selected by an essentially arbitrary 
choice. In order for AJ to become a required extreme, it must satisfy two 
conditions: 


dAJ _ 3 AJ _ q 
da 3b 


( 12 ) 
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Find a position of S(a,b) relative to the 
mode $(x,y) such that AJ becomes an 
extreme [maximum or minimum depending on 
the sign of 4>(x,y)]. 

AJ is a measure of vibroacoustic coupling. 


Figure 2-1. Illustration of a Vibroacoustic Coupling for a 

Planar Structure 
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For a general case of a mode, ^exceptions (a^ase^f very simple 

that, sofar, remains unresolved. Wit the P CO ntour of As or where Ap»As) 

r«r o? ^bW,nr, &r&Tu z rz\ e «* 

, *• A ci a h! QAtisfvina (12) may use a lengthy technique of a 

A strategy to find S(a,b) sa ^ st { 1 v "£ n i 0 parameter a is varied until an ex- 
single parameter variation; for exa p , par^ ^ second parameter b is varied 
treme where 3AJ/9a*0 is found at^a ai- The entire pr0C ess can 

and another extreme where ^/ab~0 is fou^d^ 4n l 4ccept451e convergence. 

be repeated as many times as necess y acceDtable if performed by a 

The strategy may appear faultless and, perhaps, acceptaoie it p 

computer code. 

Aside from the fact that such an ^hl^teri^ strategy were 

dynamic analysis codes, f ) W1 where conditions (12) are satisfied, may de- 
developed. The point S(ai,bi|, absolute extreme. Particularly in 

st“rt^s°Ih^re r f ^de asymmetry^ becausedby the conH gora- 

"^hSnSl Teat; any region 

containing the center of CPD. 

The problem of multiple extrenms may occu^i^the fundamental ^jf^struc- 
^.“raS can £ ^nd’even in simple ^-jn^e^'.^tKese 

sar-sssi sissi.r 

al^T: S. ££ " _ 

extreme of the AJ-coefficient for such a mode. 

An important application case occurs JStfM 

^ te be obta1ned°as fpS) two separate estimates and. perhaps, using the 
diagrams presented in appendixes A and u. 


If 

<t»(x,y) = 4»uv(u,v) = ‘buiu'frvW 
Then using form (10) for P uv (u.v) and considering dA=dudv: 

AJ = f* u (u)P u (u)du f 4> v (v)P v (v)dv 

K u **v 
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A niece-wise approximation of the Ai contour is always possible to ensure that 

“d *!-*»<*>• Then each inte9ra ' in (1 - 2) can be eva1uated sepa - 

rately: 


= j* 4> u (u)P u (u)du and AyJy = j < fc v ( v ) f M v ) dv 


(1.3) 


so that: 

AJ = A x A y J x Jy (1 * 4 ^ 

If the approximation (1.1) is possible, then the coordinates of the point 
S(a,b) where AJ becomes an extreme can be established from a separate exami- 
nation of each component function, 4> u (u) and4» v ( v )- This co nclus10n directly 
follows from (1.4) because (12) now becomes: 


3AJ = A j ^Mx) = o when d(AxJx I = 0 


8AJ , 5(A y Jy) . „ when !<Md=0 (12bl 

3lT x x 3b 3b 

When the mode shapes become sufficiently complex, the approximation by (l.D 
should be attempted, at least for the purpose of establishing the position of 
the CPD center where AJ becomes an extreme. 
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SECTION III 

SOLUTION FOR A UNIFORM CONTINUOUS BEAM WITH EQUAL SPANS 


The example presentedin .^^.^^o^ ^^^ff^ients^y^dlrect^ntegra-* 1 
one can obtain an explicit solutio numerical integration. The 

W»T5K.S -1,1 and on,, a 

few cumbersome algebraic operations are omitted. 

Variables defining a one-dimension problem are shown In figure 3-1- A " sting 
o* indefinite integrals used in the derivation is shown in figure 3-2 for 

reference. 

The origin. of the x, y 

^°TXV^r. it U ; function of Integral limits defined by the 
overlap°area Ai . The concept of the "loaded" span Is introduced as a con- 
venience for future reference and presentation of results. 


(13) 


In the x,y system: 

4>(x,y) = *P(x) = cos 2^- 

cycles of foslne functlon representlng the ™de (n - mode number for a simply 
supported beam). 

uniquely 1 — I—’ W-coefficlents^If^the 
.ft • 


''Of 


then: 


a. 


<p(x) = cos i M C 


= i. nml = ^^M 


b. 


c. 


Presented AJ-coefficients would require a multiplier <P 0 . 

Corresponding AJ/M and the response modal coordinate q=Y(AJ/M/w n ) 
woulfTcontain a multiplier W,. However the resu ing pe k 
response modal displacement (response ampTitude) Z-q««x) is indepe 
dent from * 0 , and so are all other response quantities (stresses, 

reactions, etc.). 

The above relations should be kept in mind whenever a mode is not normalized 
to the form assumed in (13). 
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A continuous beam with pinned end supports. The beam is assumed to have a 
uniformly distributed mass m (lb-secvin^) and a constant stiffness El 
(lb-in2). 


Total number of equal spans is n. 



First mode resonsance frequency (undamped): f 


^ [ IT 

21 2 V M 


(Hz) 


Origin of coordinate system i,x is located at midspan of the loaded span no. K. 


First bending mode: <P(x) = cos 3^. 

Limits on x : - ( K - l^l<x<^n - K + .1^1 


When the mode is normalized to have a maximum displacement of one unit 
(1 inch), then the generalized modal mass, M, is equal to one-half of the 
total mass: m « J. null 
2 

Define: « * A and ft * A parameters of an explicit solution. 

1 1 


Figure 3-1. Parameters for Calculation of AJ-Coefficients for a 

Continuous Beam 
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Figure 3-2. Listing of Indefinite Integrals Used in the Derivation 
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If the center of CPD is selected at point S(-a,0), then: 

u = a+x ; x = u-a ; dx = du ( 14 ) 

If the integration of (1) is performed in the local u coordinate and variables 
are transformed accordingly by (14): 

* (x,y)-<t> uv (u,v) = ^ u (u) = cos - u ~ a) 

1 


* u (u) = cos — cos + sin ffa sin ™ 

11 11 

(15) 

P(x,y) -+Puv(u>v) = P u (u) = i^l + cos 2iru ^ 

(16) 

The integrand in (1) becomes: 


*Vu)P u (u) = i. cos M j^cos EM. + cos EM cos 2 jtu J 





sin EM + sin EH cos 2wu 1 
1 1 X J 


(17) 


There are two constraints on the allowable variation of u. The first 
imposed by the extent of the structure (same as the constraint by As): 


is 



u 

T 




(18) 


The second is imposed by the extent of the 
as the constraint by Ap): 


pressure correlation length X (same 


-X < “ * _L 
21 * 1 < 21 


(19) 


Where K is the number of the span containing the CPD center 
left to right, in the direction of the +x-axis). 


(counting from 


The oveH a p between constraints (18) and (19) defines lower and 
of the integral (1), the same as the constraint provide by Ai. 


upper limits 


J ^(x.y )P(x,y )dA 

Ai 


> 

V u ) p u( u ) du = AJ 

“1 


(la) 


3-4 



KSC-DM-3265 


Where u-j is the larger (considering signs) value between: 

1^. - K + ^ and - A (20) 

u u is the smaller value between: 

l/a. + n - K + ±\ and + A (21) 

VI Z) 2 

At this point, all constituents of (1) and (la) are defined by (17), (20), and 
(21) and AJ can be calculated explicitly by using pertinent integrals from 
figure 3-2. Before this is done, there remain a few aspects that require 
clarification. These are presented next. 

Only the total value of an AJ coefficient that is both required and sufficient 
for intended application in response spectra can be calculated. Neither A nor 
J from AJ can be separated without introducing a convention of how to do it. 

In the analysis of a steady-state response to random acoustic pressures by 
means of input power spectra, a concept of Joint Acceptance, usually desig- 
nated by the symbol J 2 , is introduced. There are tables and diagrams 
depicting J2 for fundamental and higher modes of single-span structures such 
as beams, plates, etc. Calculated quantity is actually (AJ) 2 where A is 
assumed to be the total area of the structure (a single span only) so that J 2 
can be uniquely defined. In order to remain compatible with the concept of 
Joint Acceptance and in order to eliminate the dependence of J upon the span 
length 1, it is assumed that A=1 . This convention, actually a normalization, 
allows a presentation of J in the form of a nondimensional coefficient that 
has the same physical meaning as the square root of Joint Acceptance. In the 
following equations, AJ and 1J are used interchangeably to denote a product of 
two defined quantities. 

A general formula for the integral in (la) contains indeterminate terms of the 
type 0/0, when X=21 and when X=0. The last one is a trival case of zero load 
when a correspondi ng AJ=0. The case X=21 requires a different formula for J 
which does not contain indeterminate terms. 

In order to simplify an already cumbersome explicit formula for J, the follow- 
ing definitions are used: 

« -A , p-l , x - Jij- , 

c = ff( P - 2 ) t D = ff( P t Z1 (22) 

P P 
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Then, constraints on U] and u u , (20) and (21), become constraints on X and Y: 

X = MA x(-| , <* - K ; Y = MIN (| , «* + n - K +±) (23) 


where K=l,2. . .n 
For P¥=0 and 0*2 


J = ! 
1 


i)du = 


- si 


if «A U (U)P U (U] 

J u 1 

r sin ffY - sin irX + sin CY - sin CX + sin DY - sin DX | 

L 7 2C 2D J 

nira f cos irY - cos irX + cos C Y - cos CX + cos DY - cos DX | 

— L 7 2C 20 j 


COS IT g 
2 


For 0=2 


j = cos 


irY - sin jtX + sin 2*Y - sin. 2*X + 

4 2 


_ sin 

2ir 


For 0=0, J=0 


1SL j^sin 

[ cos »Y - cos irX + cos m - cos.iiX j 


(24) 


(25) 

(26) 


Explicit equations (24) through (26), definitions (22) and constraints (23) 
complete the solution of equation (1) for the case of a uniform continuous 
beam with equal spans. 

The examples in appendix A contain calculated J-coefficients ranging from the 
case of a single span beam (a trivial case included mainly for a comparison) 
to the case of a five-span continuous beam which appears to represent the 
largest number of continuous spans encountered in practice. The accompanying 
text on the plots provides an explanation of the basic plot parameters. The 
user is also referred to figure 3-1, which depicts plot parameters. 
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SECTION IV 

AN EXTENSION TO SPECIAL CASES 


The solution obtained in Section III can be extended to a few special cases: 

a. A simply supported beam with two equal cantilevers 

b. A free-free beam 

c. A single cantilever 

Extension of the existing solution requires only a few simple modifications in 
the formulation of constraints on the integral limits and an addition of a 
constant term to the equation defining the fundamental mode of case struc- 
tures. It should be noted that new solutions are approximate, although quite 
accurate, because they are based on an assumed mode shape rather than on an 
exact solution for the mode. While exact solutions exist, their formulation 
contains both trigonometric and hyperbolic functions so that the roots of 
characteristic equation which define resonance frequencies cannot be obtained 
explicitly but only by numerical methods. A check of the accuracy of an 
assumed mode solution may be based on a comparison between fundamental reso- 
nance frequencies obtained from the exact solution and an approximate solu- 
tion. Such a check shows that for a free-free beam (a particular case of a 
simply supported beam with equal cantilevers) the error is approximately 1.3 
percent. When cantilevers increase relative to the midspan, in the limiting 
case of the midspan length approaching zero (a case of two equal cantilevers 
with a common fixed end), the error increases to a maximum of 4.1 percent. 
Thus, within the entire range of the approximate solution, the error from the 
approximation remains within acceptable limits. 

Figure 4-1 shows a new configuration and the variables defining the extension 
of the previous solution to special cases considered in this report. The 
assumed fundamental bending mode is: 

4>(x) = cos *21 - sin*7 =4 >sb(x) - F 

where 

4> sb (x) = cos 121 


is the first bending mode of a simply sup- 
ported beam with the span L. The form of 
4>SB(x) is identical to that of v>(x) in (13) 
when 1 is substituted by the new variable 
L=l+2d = the total length of the beam (center 
span and two cantilevers). The solution for 
J=JSB» when the mode is 4>s B (x), is 
already provided by equations (22) through 
(26) after the substitution of 1 by L is made. 
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A simply supported beam with two equal cantilevers, beam is assumed to have a 
uniformly distributed mass M (lb-sec 2 /in 2 ) and a constant stiffness El 
(lb-in 2 ) 



First mode undamped resonance frequency (by Rayleigh energy method) is: 

f * -JL_ IE |"l - JL s 1 n »7 + Zsin^iryl = .* [EL (Hz) 

2L 2 V m L » J 2L 2 v 2M 

Generalized modal mass: M = - 1 sin *7 + 2 sin 2 * 7 j 

Define: « = A and 0 = 


4-2 


Figure 4-1. Configuration and Parameters for Calculation of 
AJ-Coefficients in a Case of a Simply Supported Beam With Two 

Equal Cantilevers 



KSC-DM-3265 


is a constant in the case of a particular 
structure. The new variable, 7, defines a 
ratio between the cantilever length, d, and 
the total length, L. As noted in figure 4-1, 
each value of 7 defines a different struc- 
ture. Consequently, J-coefficients for dif- 
ferent structures are not directly comparable. 


4>(x) in (27) represents a superposition of two inodes: $sb(*) < a bending 
mode) and F (a translation mode). Both modes have the same resonance fre- 
quency. Consequently, the J-coefficient for *(x) becomes a sum (a super- 
position) of two coefficients, J$b a "d Jf, that correspond to ^SBU) an J 
F modes. Such a summation is valid, provided that Jp is referenced to the 
same A=L as J$b is referenced to in equations (24) and (25): 


j = J SB - Jp 

The same result as (28) can be obtained by using <t(x) 
the same procedure outlined in equations (14) through 
is simplier. 


(28) 

in (27) and by repeating 
(26). The chosen route 


This paragraph provides an explicit definition for the generalized modal mass, 
M for the mode 4>(x). The component mode 4 >sb(*) is normalized to the unit 
maximum displacement at x=0. The amplitude of the second component mode, F, 
is a function of parameter 7. Unlike the case of a continuous beam with equal 
SD ans where M is equal to one-half the total mass, the generalized modal mass 
and the resonance frequency of a cantilevered beam become functions of the 

parameter y: 


r 1/2 

M = I 4*2(x)mdx = ££ |"l 
J -L/2 


£ sin»7 

IT 


+ 2 sin 2 ir7 




(29) 


If an "equivalent" uniformly distrubuted mass 


M 



is introduced, then the fundamental resonance frequency of a cantilevered beam 
attains the same form as the fundamental resonance of a simply supported beam: 



New parameters in (22) are now defined in terms of L: 



4-3 


KSC-DM-3265 


The form of parameters C and D in (22) remains unchanged. The constraints on 
X and Y become: 


X . MAX (-f ,« - 1), Y - MIN(|, « + 1) 


(23a) 


In terms of the parameters of (22a) and constraints of (23a), equations (24) 
through (26) now define J$g. The equation for the remaining term Jp in 
(28) is: 


,“u 


J P u (u)du = siniry f~ Y-X + _8_ 

( sin hi - 

sin hi 

J U] L 2 4ir 


0 )\ 


(31) 


Equations (28), (22a), (23a), (24) through (26), and (31) complete the ex- 
plicit solution for J-coefficients of a cantilevered beam and a free-free 
(7=0.224) beam. A comparison between J-coefficients for cantilevered beams 
with the same total length L is not meaningful because the values of corre- 
sponding M s in (29) and f's in (30) are variable functions of the parameter 
7- 


The solution for a single cantilever with the span d, for example the left 
cantilever in figure 4-1, is obtained from the proceeding solution by setting 
7=1/2 and by constraining the Y-extent of the CPD load to: 



The corresponding generalized modal mass and the fundamental resonance fre- 
quency are (approximate error 4.1 percent in f) from (29): 

M = M [3 - i] ** 0.227md 


The exact solution yields M=0.25md (error » 9.3 percent) from (30): 



(The exact coefficient is 0.56 instead of 0.583.) Then, using (28) and L=2d: 

AJ = LJ = L(J sb - J F ) = 2d(JsB - Jf) = dJ c (32) 

and J c = 2 (Jsb - Jf) (28a) 

Here, J c is the J-coefficient for a single cantilever, referenced to the 
cantilever span d. All parameters in (22a) entering computations of Jcr and 
Jp) remain defined by L=2d. Thus, horizontal plot axis, 0, defines a ratio 
between the pressure correlation length and two cantilever lengths. These 
special cases are shown in appendix B. 
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SECTION V 

SUMMARY AND RECOMMENDATIONS 


The existing dynamic analysis codes do not have the capability to utilize re- 
cently developed methods involving response spectra and the distributions of 
correlated pressures in the analysis of structural responses to random acous- 
tic loads. The implementation of such a capability within an existing code 
would require a substantial modification of the code that may be far beyond 
what the code structure allows to be accomplished by a user. Perhaps, only 
the originator of the code may have the capability to implement such modifica- 
tions. 

A different route should be considered. This route would use either an exist- 
ing FEM code or a modal analysis test to establish the modal parameters of a 
structure. The capability to have a dual source of modal parameters is essen- 
tial for a test/analysis correlation and for a dynamic analysis of structures 
where only the test can establish modal parameters. A separate code should be 
developed for the analysis of structural response to random acoustic loads 
occurring during a launch. The code would input a structure geometry, modal 
parameters, pressure correlation lengths, and response spectra. The input of 
geometry and modal parameters should be made by the transfer, via a magnetic 
tape or other media, from the modal analysis code to the new response analysis 
code. For a general case of a three-dimensional structure, a corresponding 
three-dimensional definition of the contours of correlated pressure distribu- 
tions over the structure surface must be obtained from the ellipsoid of pres- 
sure correlation lengths. The original grid of nodal points, which was gener- 
ated by the requirements of modal analysis, may have to be either replaced or 
appended by a high-density grid where the density is governed by the non- 
uniformity of the correlated pressure distribution and by the requirement to 
have a reasonably small delta-A associated with each grid point to ensure the 
accuracy of a numerical integration. Modal vectors will also require a trans- 
formation from the modal analysis grid to the high-density pressure distribu- 
tion grid. Then, the estimate of AJ extremes and peak response modal coor- 
dinates can be performed in the high-density grid coordinates. Computed modal 
coordinates would be either input back into the modal analysis code for the 
calculation of peak responses (stresses, displacements, etc.) or the same task 
may be accomplished within the response analysis code, in which case addi- 
tional inputs, such as matrices of internal stresses due to normal modes, must 
be obtained from the modal analysis code. 

In the above outline, the optimization of the load position leading to an 
absolute extreme of AJ-coefficients and modal coordinates may never become 
entirely automated, and the problem is not expected to be solved for a general 
case. Because each mode requires a different load position for the extreme 
response in that mode, the initial positioning of the CPD center may have to 
be left to the judgement of the user, followed by some limited programmed 
optimization. A graphic display of the mode and an interactive input of the 
initial load position would certainly facilitate the choice of input. 
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In this scheme, a separate response analysis code provides a definite advan- 
tage before a general modal analysis code. However, the communication prob- 
lems between the two codes involving data transfers, stops and restarts, are 
not expected to be resolved without the involvement of personnel who origi- 
nated and presently maintain the general modal analysis code. 

A development of the capability to perform a response analysis for a general 
case of a structural configuration would require substantial effort and time. 
Perhaps, it should be accomplished in steps of a gradually increasing com- 
plexity by limiting the initial development to only certain types of struc- 
tures, such as planar structures, etc. 


Meanwhile, in order to satisfy the immediate needs of design groups, designers 
must resort to approximations, such as breaking down a structure into simpler 
components for which adequate solutions either already exist or «n be ob- 
tained by means of simple programs and diagrams like those presented in this 


report. 
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APPENDIX A 

J -COEFFICIENTS FOR ONE, TWO, THREE, FOUR, AND FIVE SPAN BEAMS 
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APPENDIX B 

J-COEFFICIENTS FOR A CANTILEVERED BEAM, A FREE-FREE BEAM, 
AND A SINGLE CANTILEVER 
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B-6 


BETR=CORLEN/ ( FREE BERM LENGTH)x 10 

Horizontal axis , BETA , defines a RATIO between PRESSURE CORRELATION LENGTH and the TOTAL BEAM LENGTH. 
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APPENDIX C 

EXAMPLE PROBLEM FOR A THREE-SPAN CONTINUOUS BEAM 


This appendix contains numerical examples illustrating the application of the 
presented theory to a computation of response of a three-span continuous beam. 

In the presented examples, the response spectra, pressure correlation lengths, 
and modal parameters are assumed to be known. Under this assumption, the 
response computations are shown. Most of the text is devoted to step-by-step 
explanations of the employed procedure. 
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The example of a three-span continuous beam is shown in figure C-l. 


Beam parameters are: 

Number of spans 
Span length 
Beam width 
Uniform mass 

Material modulus of elasticity 
Moment of inertia 
Section modulus 
C.l FIRST SINUSODIAL MODE 
The resonance frequency (undamped) is: 

fi = _JL_ */H = 30 Hz 
1 21 2 

u>l = 2irfi = 188.5 rad/sec 


n = 3 

1 = 58.0 in (4.833 ft) 
b = 7.2 in 

m = 1.6814 x 10-3 lb-sec 2 /in 2 
E = 29.0 x 10 6 lb/in 2 
I = 0.2393 in 4 
S = 0.282 in 3 


co ^ = 3. 553 x 10 4 (rad/sec) 2 

The generalized modal mass (mode normalized to the maximum unit displacement 
at the centerline span) is: 


M = i. nml = 0.1463 lb-sec 2 /in 
2 

The area exposed to acoustic pressures (for application in the diagrams of 
J-coefficients) is: 

A = bl = 417.6 in 2 

The response spectrum ordinate for an assumed 2-percent damping (see figure 
C-2) at f i=30 Hz is: 

2 

Y = q/(AJ/M/co^) = 6.5 
This is a maximum value from 12 launches. 

Note that the median value of Y is around 4. The assumed value, Y=6.5, is the 
maximum design value. 
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Figure C-l. Test Model Data 
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Figure C-2. Lift-Off Response Spectra to Pressure on the Fixed Service 
Structure (FSS) at 155-Foot Elevation (Sensor KSAPB002A) 
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The pressure correlation length, X, depends on the direction of the beam (see 
figure C-3), whether the beam is vertical or horizontal. 

C.1.1 CASE OF A VERTICAL BEAM. 

At fl = 30 Hz, X v « 19 ft 

0 V = = 3.93 

From the diagrams of J-coefficients for a three-span continuous beam, maximum 
J v at 0 V occurs when the CPD center is locted in span no. 1: 

maximum J v = 0.43 (see figure C-4) 

Note that if the CPD center were located in span no. 2, the error in the 
response estimate would exceed a factor of 4 (see figure C-5). 

2 

AJy/M/ojj = 417.6 x 0.43/0.1463/3.553 x 10* = 3.455 x 10-2 
The response modal coordinate (modal participation factor) is: 

q lv = Y ^AJ v /M/a^ j = 6.5 x 3.455 x 10~ 2 = 0.2245 

Because the mode was normalized to the maximum unit displacement, 

<117=0.2245 inches is the maximum vibration amplitude at the centerline span. 

C.1.2 CASE OF A HORIZONTAL BEAM. 

At fi = 30 Hz, X h = 75 ft 

J3 h = LH = 15.5 

Because of a large 0h» the location of the CPD center yielding a maximum Jh is 
not critical, whether it is in span no. 1 or no. 2. 

For the CPD center in span no. 1, = 0.563 (see figure C-6). 

For the CPD center in span no. 2, = 0.58 (maximum) (see figure C-7). 

AJh/M/tOj = 417.6 x 0.58/0.1463/3.553 x 10 4 = 4.66 x 10-2 

q lh = Y AJ h /M/c^l = 6.5 x 4.66 x 10-2 = 0.3029 
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Figure C-3. Lift-Off Acoustic Pressure Correlation Lengths on FSS 
Between 125- and 155-Foot Elevation 
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Figure C-4. J-Coeff icients for a Three-Span Continuous Beam, Correlated 
Pressure Distribution Center in Span No. 1 
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Figure C-5. J-Coefficients for a Three-Span Continuous Beam, Correlated 
Pressure Distribution Center in Span No. 2 
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Figure C-6. J-Coefficients for One-, Two-, Three-, Four-, and Five-Span 
Beams, Correlated Pressure Distribution Center in Span No. 1 
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Figure C-7. J-Coefficients for Two-, Three-, Four-, and Five-Span Beams 
Correlated Pressure Distribution Center in Span No. 1 
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C.1.3 FIRST BENDING MODE SUMMARY. Modal coordinates, either 9i v or are 
the multipliers of the normal mode no. 1 stress matrix. For the first mode, 
^=sin(irx/l ) , stresses can be computed explicitly without resorting to the 
FEM. Thus, the stress matrix bending moment is: 

m ■•«{£■« ( 

2 

M 51 (at x = 1/2) = El = 2.036 x 10 4 lb-in/ (normal mode amplitude) 

For a vertical beam, the maximum response bending moment at midspan 
(x=l/ 2 ) is: K 

5M r blv = q lv M bl = 0*2245 x 2.036 x 10 4 = 4.571 x 103 in-lb 
and the maximum unit stress is: 

5 rblv = Mrblv/ S = 4 * 57 1 x 103/0.282 = 16,209 psi 
For a horizontal beam: 


sin (ttx/1 ) 


5 rblh = q lh Mbi/S = 21,870 psi 
C.2 SECOND BENDING MODE (ANTISYMMETRIC) 


An approximate estimate of beam response in this mode is possible using 
J-coefficient diagrams presented in this report. Note that in the following 
example, J-coefficients are not those of a three-span beam. 9 


Modal parameters are taken from the FEM [NASA Structural Analysis (NASTRAN) 1 
output: 


Resonance frequency f 2 = 38.45 Hz 

Generalized modal mass M 2 = 9.599 x 10~2 lb-sec2/in 

(mode normalized to maximum unit displacement in spans no. 1 and no. 3 ) 

gj ^ = ( 2irf 2 ) ^ = 5.836 x 10 4 
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C.2.1 CASE OF A VERTICAL BEAM. 

At fz = 38.45 Hz, X v 2 » 15 ft 


0 v2 = hi = 1§- = 3.10 

1 4.833 

Locate the CPD center at midspan of span no. 1 and draw (in scale) the 
correlated pressure distribution, 


P(x) = I 
2 


1 + cos 

X V 2 


, considering that 


= ^v2 i = i. 55 i as shown in figure C-8. 
2 2 


Write the expression for AJ, paying attention to integral limits: 

- X v 2/2 /• 1/2 /*^ v 2/2 

AJ = I V» 2 (x)P(x)dx = I V> 2 ( x )P(x)dx +1 * 2 U)P(x)dx 
^- 1/2 - 1/2 J 1/2 

Note that: 

a. For 1/2 < x< X v2 /2, ^(x) is small compared to ^(x) in span no. 1. 


r X v2/2 

b. I V> 2 (x)P(x)dx will have the sign opposite to 
J 1/2 

r 

I ^ 2 (x)P(x)dx 
- 1/2 
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Therefore, a conservative estimate of AJ can be obtained by calculating only 
the first of the two integrals: 

r 

AJ * I 

- 1/2 

Note that in span no. 1, is similar to the first mode of a simply 

supported beam with span 1. Therefore, one can compute AJ and J from the 

diagram (see figure C-9) for a simply supported beam at 0 V 2=3.1O , J V 2~0.607. 

Corresponding A=A V 2 is the same as for the first mode. 

AJ v 2/M2/<4 = 417 - 6 x 0-607/9.599 x 10-2/5.836 x 104 = 0.0452 

The response spectrum ordinate for an assumed 2-percent damping (same plot as 

for the first mode) at 38.45 Hz is Y2=9.5. Y?=9.5 is a 12-launch maximum. 

The next highest Y is around 6.5, approximately 32 percent lower. 

92v = Y2 AJV 2 /M 2 /W 2 ] = 9 : 5 x 0.0452 = 0.430 

Response in the second mode exceeds that in the first mode. This rather 
unusual case is a consequence of small high Y 2 at launch STS-7. 

From the modal stress matrix (NASTRAN output) the bending moment at the 
midspan of span no. 1 is: 

Mb2 = 2.385 x 104 lb-in/ (normal mode amplitude) 

The response bending moment is: 

M rb 2v = q 2v M b2 = °.430 x 2 * 385 x 104 = 1,025 x 104 in_1b 
The maximum unit stress is: 

6 rb2v = = 36,356 psi 

C.2.2 CASE OF A HORIZONTAL BEAM. 

At f2 = 38.45 Hz, X h 2 * 67 ft 

/J h2 = = 13.87 
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Figure C-9. J-Coefficients for a Single-Span Simply Supported Beam, 
Correlated Pressure Distribution Center Within the Loaded Span 
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If the same steps are followed as for the vertical beam and the CPD is drawn 
relative to the mode, considering that: 

^2 = 6.941, 

one finds that the entire asymmetric mode is "loaded" by a nearly uniform 
distribution. If the distribution were exactly uniform, the mode could not 
be excited since 


f 

J ^(x)dx = 0 

An approximate estimate of AJ can be made by a direct computation. The CPD 
center is positioned at the left support of span no. 1 because this position 
results in the most nonuniform load and maximum AJ. This is shown in figure 
C-10. 

At x = 1/2, |l* = iL = -IL. = 0.2265 
*h2 X h2 ^h2 

P(l/2) = I Cl + cos (0.2265)]= 0.9872 = Pi 


At x = 51/2, fSi = |I_ = 1.1325 
x h2 *h2 

P (§l\ = I [1 + cos (1.1325)] = 0.7122 = ?2 
\2 / 2 

It is assumed that: 

r 

J v>(x)P(x)dx = U 

The assumption is conservative since the actual value of this integral is 
negative and it would, if computed, decrease the total estimate of AJ. 

r 31 

(AJ)h = I bv>(x)P(x)dx KP 1 -P 2 ) 


= 7.2 x 1 x 58 x (0.9872 - 0.7122) = 76.6 in 2 
3 


( AJ )h/M2/cc>2 = 76 - 6 / 9 - 599 x 10" 2 /5.836 x 10 4 = 1.367 x 10“ 2 
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Figure C-10. Second Bending Mode in the Case of a Horizontal Beam 
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^2h = ^2 ^(AJ)h/M2/w2 


= 9.5 x 1.367 x 10-2 = 0.1298 
The response bending moment is: 

M r b2h ■ q 2h M b2 = 0-1298 x 2.385 x 10* = 3.097 x 103 in-lb 
The maximum unit stress is: 


6 r52h = M2h = 10,981 psi 

C.2.3 SECOND BENDING MODE SUMMARY. Estimates of modal coordinates are 
approximate and conservative. For a vertical beam, J was computed from the 
diagram (see figure C-9) for a single span beam. Such approximation was 
possible because is small and the contribution to AJ from the center span 
is negligible [small For a horizontal beam, AJ was estimated directly 

by calculating the integral definition of AJ. In this computation, the fol- 
lowing general theorem was used (see the following illustration). The theorem 
states that: 


* 

J 


^(xjPitxJdx = A^jPiUcg) 


The integral is equal to the area, Ay>, under ^(x) times the ordinate of the 
linear function, Pi(x C g), taken at the abscissa, x C g, corresponding to the 
location of the center of gravity of A^. 
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function in the interval 


Any function having the same 
sign in the interval o<x<l 
and "center of gravity" at 

x cg- 


Area under ^(x) 



(x)dx 


Position of the 
"area center" 


x cg 



v>i(x)xdx 


Proof of theorem: 


Pl(x) = Pi(o) 




Pl(x)<Pi(x)dx = P( 



1 - 


^(xldx 



1 


^(x)dx 

o 



^l(x)xdx 
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« Pl(o) [A^ - A^-^S.] " A/i(o) [l - 


= A^pP ]_ ( g ) 

C.3 THIRD BENDING MODE (SYMMETRIC) 

This mode is a result of a mathematical solution of the eigenvalue Problem. 

J-n 

check the program) can be computed (approximately) using the following 
procedure: 

a. Locate the CPD center at the center of span no. 2. 

b. Use a procedure similar to that in C.2.2, applying the theorem in 
C.2.3. 

c. When calculating the area of modal displacements in spans no. land 
no. 3, take into consideration modal amplitudes in these spans (t y 

are =£ 1.0). 

d. Within each span, the correlated pressure distribution function 


P(x) *|[ 


1 + cos 


2wx 
X J 


shall be substituted by a linear approximation so that the theorem in 
C.2.3 will apply. 


C.4 FOURTH BENDING MODE (ANTISYMMETRIC) 


This mode shape is the same as the first bending mode of a : si: an ^ 
the snan l/i=l/2=29 in=2.42 ft. The mode is shown in figure C-12. Because tne 
mode 5 ? ^anti symmetric, the critics! position of the CPD center is elwevs in 
the first snan As 1 onq as X/2 (X ■ either vertical or horizontal PCL) does 
not exceed 5xl 4 =12.08 ft, existing diagrams of J-coefficients derived for a 
five-span beam can be used. 
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Figure C-12. Fourth Bending Mode 
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In this case, 

A 4 = bl 4 = bl/2 = A/2 = 208.8 in 2 


Mode no. 4 parameters (from NASTRAN output) are: 

Resonance frequency f 4 = 120 Hz 

Generalized modal mass M 4 = 0.1617 

w 2 = ( 2irf 4 ) 2 = 5.685 x 10 5 

Because the mode resonance, f 4 , exceeds the frequency range of available PCL's 
and response spectra, the required values of X's and Y 4 will be extrapolated/ 
assumed. From a consideration of available power spectral densities (PSD's), 
the value of the response spectra ordinate at 120 Hz is Y 4 =12 (±10 percent). 

C.4.1 CASE OF A VERTICAL BEAM. 

Extrapolated X v4 = 4.8 ft 

1 4 

From the diagram of J-coefficient for a five-span beam (see figure C-13), at 
0 V 4=2, J v 4=0-538. 

A 4 J v 4 /M 4 /^ = 208.8 x 0.538/0.1617/5.685 x 10 5 = 1.222 x 10 ‘ 3 

21 

Q 4v = Y 4 A 4 J V 4 /M 4 /to 4 = 12 x 1.222 x 10"3 = 0.0147 

From NASTRAN output, Mb 4 =-8.416xl0 4 at approximately one-quarter of the span 
("element 4" in NASTRAN ID). The bending moment at the midspan in this mode 
is zero. The response bending moment at 1/4 is: 

M rb4v = 8 * 416 x 104 x °-°l 47 = l * 234 i 0 - 115 
The maximum unit stress is: 

J rb4v = 5 !eMv = 4,376 psi 


C-24 



KSC-DM-3265 


c 

rO 

O. 

<A 


O' 

U 

5 i 

o to, 

^ { i 

w or I 
Pi in! 

• Si 


m 

CD 

a 

o 


2T SI 
»-» ■ 

2 

R m! 

0 cm. 

«— • i 

ii: 

cn 

1 ■ i in 
> ~ 
»-» R 
Li- cr 


in 

IN. 

03 

I 

(C 


CO 


c 

<D 


a 


Ai 


CD 

o 

o 


a: 

x 

a. 



CSJ 

C3 

X 

:z 

a: 

CL 

cn 

N 

UJ 

_J 

(Y 

O 

o 


II 

a: 

i — • 
UJ 
CD 


( 


aouB^daooy ^uiof 


}ooj - bs) 


f 


* 


TO 

TO 


I 

B 


r 

t 

s: 

P 


TJ 

4 

0 

3 

3 

* 


C-25 


Figure C-13. J-Coefficients for a Five-Span Continuous Beam, Correlated 
Pressure Distribution Center in Span No. 1 
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C.4.2 CASE OF A HORIZONTAL BEAM. 

Extrapolated = 15 ft 

0 h4 .i44. 6.2 

u 

From the diagram of a five-span beam (see figure C-13) with the CPD center in 
span no. 1, = 0.32 


A 4 J ^ 4 ^ 4 / 0 )^ - 7.27 x 10 -4 


q 4h = y 4 


A4Jh4/M4/oj 


= 8.72 x 10-3 


The response bending moment at 1/4 is: 

M rb4h = q 4h M b4 = 734 1n “ lb 
The maximum unit stress is: 

5 rb4h = 2603 P S1 ' 


C.4.3 FOURTH BENDING MODE SUMMARY. Response was computed using 
J-coefficients for a five-span beam. The response is diminishing. No higher 
modes need to be considered. 


C-26 



